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Abstract 

A large fraction of SN1987A electron antineutrino events has been recorded in the first 
second. We study how this observation fits into the conventional paradigm for neutrino 
emission, and show that there is a 3. 2a hint for an initial accretion phase. This phase 
involves a large fraction of the energy emitted in neutrinos and antineutrinos, about 20% 
or larger. The occurrence of neutrino oscillations strengthens these inferences. We discuss 
why three flavor oscillations with normal mass hierarchy are completely acceptable, whereas 
oscillations with inverted mass hierarchy require more troublesome interpretations, if #13 
is above 0.5 — 1°. 
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1 Context and motivation 

The first second after a gravitational collapse is a moment of crucial importance. The very intense 
initial neutrino luminosity (denoted as 'accretion' in the following) is expected to have a non-thermal 
character [I], and it is thought to be the key to understand the subsequent explosion of the star [21 [3] . 
In this connection, it is very interesting to note that about 40% of the SN1987A events have been 
recorded in the first second: 6 out of 16 in Kamiokande-II [3J, 3 out of 8 in 1MB [5] and 2 out of 5 
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in Baksan [6] (the comparison is made with the number of events recorded in a window of T = 30 
seconds). This theoretical expectation and this experimental fact motivate us to analyze quantitatively 
the first second of SN1987A. 

In our calculations we will largely follow Lamb and Loredo [7], who included for the first time a 
description of the background and of the time-energy distribution of the events. We will point out 
in the Appendix the technical points where we depart from their analysis. We will discuss the role 
of neutrino oscillations in the interpretation of the first second of SN1987A neutrinos and stress their 
importance. 



2 Formulation of the problem 



The parameter that we aim to study is the fraction of energy that was emitted in the non-thermal 
phase of neutrino emission, that occurred in (a fraction of) the first second. In formulae: 

fa = ~~p. p (1) 

where the suffix a stays for 'accretion' (or non-thermal phase) and the suffix c for 'cooling' (or thermal 
phase). Strictly speaking, this fraction cannot be reconstructed completely from the observations, 
since-in a very reasonable approximation-we saw only electron antineutrinos. Thus, in order to fulfill 
the task, we must rely on some theoretical assumption here. We will assume unless stated otherwise 
that 

( £ a = 2 x £ a (v e ) ^ 
\ £ c = 6 x £ c {v e ) 

namely, we declare that the ratios between the total energy radiated £, and the energy radiated in 
electron antineutrinos £{v e ) is 2 for the accretion and 6 for the cooling phase. These numbers are 
crucial for the interpretation of the result, so let us pause to discuss them before continuing. 1) The 
first fraction describes the assumption that during accretion only v e and v e are radiated in equal amount 
due to ep — > nv e and e + n — > pv e . It does not seem implausible that the v e are even more abundant than 
the v e during accretion, or that some other species of (anti)neutrinos are also radiated. In other words, 
the factor 2 could be an underestimation, and thus the fraction f a that we estimate from the electron 
antineutrino data can be regarded as a reasonable value (or lower bound). 2) The second fraction is 
just the usual and often adopted "equipartition" hypothesis. Only if there is a large amount of energy 
radiated in non-electronic neutrinos (say, £ c /£ c {i> e ) ~ 10 or larger) it could be possible to diminish 
significantly the ratio f a and modify somewhat the conclusions that we will describe later in this paper. 

In the rest of this section we discuss the tools that we use for a quantitative evaluation of £ a {v e ) 
and £ c (pe) from SN1987A neutrino data: first of all, we give a description of the antineutrino flux, then 
we model the expected signal rate, and finally, we discuss the likelihood function that we adopt. 

2.1 Parameterized antineutrino flux 

Let us describe the adopted form of the parameterized neutrino fluxes (differential in the energy) . We 
follow the one proposed in [7]: 



1 7TC 
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where g(E, T) = E 2 / (l+exp(£'/T)). This describes an isotropic emission from a distance of D = 50 kpc. 
The first term is given by the product of the number of targets (neutrons, with Y n = 0.6) in the accreting 
mass M a , times the thermal distribution of positrons g (with average temperature T a ), times the cross 
section of positron interactions, that increases quadratically with E v and thus gives a non-thermal 
character to the emitted v e . The second term is instead a standard black body emission from a sphere 
with radius R c . The time scales of the two processes of accretion and cooling (r a and r c ) appear in the 
functions: 

£ W - i+i/ (0 .5 S) (4) 

T(t) = T c exp[-t/(4r c )] 

Thus, we have 6 parameters (3 for each phase): M a , T a and r a for accretion; R c , T c and r c for cooling. 
For any set of values of these parameters, it is straightforward to calculate the energy carried by 
antineutrinos during accretion and during cooling, denoted by £ a {ve) and £ c (^e)0 In this way, using 
eq. El we can evaluate the value of f a . 

2.2 Signal rate 

The signal rate, differential in time, positron energy E e and cosine of the angle 9 between the antineu- 
trino and the positron directions is: 



S{t,E e ,cos9) = N p -^—(E v ,cos9) ? ]d (E e ) &(cos0) *B B (t,E u ) ^ (5) 
a cos 9 dE P 



where N p is the number of targets (free protons) in the detector, a is the v e + p — > n + e + (inverse 
beta decay) cross section, % the—detector dependent-average detection efficiency, £<j is the angular 
bias =1 for Kamiokande-II and Baksan whereas for 1MB £d(cos#) = 1 + 0.1 cos 9 [8], finally <E>p e (the 
electron antineutrino flux differential in the antineutrino energy E v ) is as in eq. O The expression of 
the antineutrino energy E v as a function of E e and cos 9 is given in the Appendix. Later we will use 
the shorthands by S(t,E e ) = J S(t, E e , cos 9) dcos9 and S(t) = J S(t,E e ) dE e . 

2.3 The assumed likelihood 

We estimate the parameters by evaluating: 

X 2 = -2 J2 log(£ d ) (6) 

d=k,i,b 

where is the likelihood of any detector (k,i,b are shorthands for Kamiokande-II, 1MB, Baksan). 
We use Poisson statistics. Dropping constant (irrelevant) factors, the unbinned likelihood of the three 
detectors are: 

C d = e - fdI -^ {t+td)dt ^e s ^^ Jsit^t^E^G^dE} (7) 

i=l 

where of course the dependence on the 6 model parameters is contained in S. Each detector saw JV^ 
events; their time, energy and cosine with supernova direction are called U, Ei and c, (i = l...Nd). 



1- We get £ a = 4.14M a T*T a <p and £ c = 3.39 IQ- a R 2 c T*t c measuring £ a , c in foe (= 10 51 erg), M a in M , R c in 
km, T a ^ c in MeV and r ajC in seconds; ip = dx exp(— x 10 )/(l + x r o /0.5) ~ 0.6 in the relevant r a range. 



3 



The time is counted from the first detected event; namely, we set t\ = for all detectors. The integral 
over the time in the first exponential factor is performed from the moment when the first neutrino 
reaches the Earth t = —tj (where td > 0), till the end of data taking, t = T with T = 30 s. The 
values of the 3 new parameters td, called 'offset times', are estimated together with model parameters 
by fitting the data (and, a posteriori found to be small) since the measurement of the absolute times 
in Kamiokande-II and Baksan are not reliable. In 1MB, the live time fraction is fd = 0.9055 and the 
dead time is = 0.035 s, whereas for the other detectors fd = 1 and Td = 0. The specific background 
rate is Bi = B(Ej) as discussed in [9] and in the Appendix (we denote by B(E e )/2 the background 
distribution, differential in time, energy and cosine-the factor 1/2 is for an uniform cosine distribution). 
The Gaussian distribution Gi includes the estimated values of the energy Ei and the error of the energy 
5Ei for any individual event; the inclusion of the error on the measurement of cos 9 does not add 
significant information, and the relative time of any event is precisely measured. 

3 Results of the analysis 

In this section we present the results of our statistical analysis. For reasons of clarity and for a more 
direct comparison with previous results, we will ignore the occurrence of neutrino oscillations in the 
first part of this section; in practice, the results that we obtain can be regarded as the best 'effective 
antineutrino flux' that describes the data. In the second part of this section, we will consider the 
modifications due to the occurrence of neutrino oscillation, and discuss why certain cases with inverted 
mass hierarchy are disfavored by SN1987A data. 

3.1 Fraction of energy in the accretion phase 

First of all, we searched the global minimum of the x 2 ■ The best fit of the initial accreting mass 
M a ~ 5 Mq is not physically plausible, since we expect that the mass exposed to positrons will be at 
most the mass of the whole outer core ~ O.6M . Thus, we analyzed the x 2 as a function of M a . For 
a wide range of values of M a we marginalized away the remaining 8 parameters by minimizing the x 2 
(normalized as usual to the best fit point) and obtained in this way fig. [TJ. From fig. [T] we conclude that: 

(i) Although larger masses would fit the data better, a completely reasonable value of the outer core 
mass, M a = 0.5 Mq, is not significantly disfavored in comparison to the best fit value. This justifies 
the fact that we adopt this value for reference in a large part of this paper, rather than the best fit 
value H 

(ii) There is a significant hint for an accretion phase. In fact, let us ask whether the improvement going 
from M a = (no accretion) to M a = 0.5 Mq is simply due to the presence of the two new variables 
T a and r a (i.e., let us perform a likelihood ratio test with two degrees of freedom). We find that we 
can reject the null hypothesis in favor of the hypothesis that accretion occurred with a significance of 
a = 1.2 x 10~ 3 , i.e., 3.2a in Gaussian language. 

(Hi) Even small values of M a are able to improve significantly the fit to SN1987A data. In fact, when 
M a decreases, the best fit value of the temperature T a increases, as one can read clearly from table HJ 
Keeping in mind eq. O one understands that the number of events during accretion scales as M a x T®~ s , 
so that small changes of T a are able to 'compensate' the decrease in M a . A similar argument applies 
also to the values of the accretion energy S a and of the energy fraction f a , that, as we can see from 

2 This argument and this conclusion is in agreement with what found in [7]- 
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Figure 1: Values of x 2 a function of the initial accreting mass M a in units of solar mass. Also 
indicated the corresponding values of f a for the reference case M a = O.5M and for the cases 
selected in tabled For M a = (no accretion) x 2 — 14.7. 

table Q] and fig- E react slowly to important changes of M a . 

(iv) Finally, a value smaller than M a = O.OIMq is disfavored at about 99% in comparison to the one 
we selected. 

Now we consider the outcome of the fit in the point M a = 0.5 M Q . The best fit of the offset times 
is zero. Their 1 sided, la errors (obtained by integrating the marginalized likelihood) are: 



similar results remain valid also with oscillations. The values of the 6 astrophysical parameters that we 
find from our statistical analysis along with the la errors (obtained instead by a conventional, A% 2 = 1, 
Gaussian procedure) are: 



The temperature of the electrons T a is pretty low; this outcome of the fit is simply due to the fact 
that the early Kamiokande-II events have low energy [9] . The duration of the accretion phase is rather 
close to the expectations, about half a second, [21 E] and [ID]. Coming to the parameters of the cooling 
phase, we see that the radius of neutrino-sphere is rather similar to the radius of the neutron star, as 
expected (see e.g., [H]). The temperature T c implies an initial average energy 3.15T C and an average 



Ai K II = 0.09 S, AtlMB = 0.31 S, AtBaksan = 0.25 s 



(8) 
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30 
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1.6 
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12 


5.5 


4.1 


36 



Table 1: The astrophysical parameters of neutrino emission defined in eqs. [3] and [7} calculated 
for selected values of M a . In boldface the derived quantities defined in eqs. U\ and® 



value 3/4 lower, namely, 12.9 ± 1.9 MeV, which compares well with the expectations [121 113| . Finally, 
the duration of the cooling phase is brief, e.g., when compared with the 20 seconds estimated in [1] 
though this leads to a reasonable value of the total emitted energy 2.4 x 10 53 erg, that compares well 
with the one expected for neutron star formation |14jl 3 1 

With the parameters in eq. [9] it is straightforward to calculate the amount of energy emitted during 
accretion, during cooling, and the energy fraction f a : 

{ £ a = 6.1 x 10 52 erg . oy . . 

\ c 10 m53 =^ /a (no OSC.) = 26% (10) 

£ c = 1.8 x lO 53 erg v ' v ' 

It is interesting to note that the central value of f a estimated from the fit is two-three times larger than 
the expected one [15] (see also [16]). Other values of f a are considered in table [Hand correspond to the 
points marked in figure [TJ 



3.2 Impact of neutrino oscillations 

We will assume that the temperature of the muon and tau antineutrinos (that are implied in the cooling 
phase) are in a fixed ratio with the v e temperature. Following [T5] , we will take as default value 

T(p r )/T{9 e ) = T{9„)/T{p e ) = 1.2 (11) 

As stated above, we also assume that in the accretion phase muon or tau (anti) neutrinos are absent 
or very rare and discuss the role of this assumption later. 



Formalism Considering oscillations among the usual three neutrinos, the observed electron antineu- 
trino fluxes are: 

= P^ e + {l- P) ^ (12) 

3 The total number of events seen / expected in 30 seconds (in brackets the background events) is: 
N KU = 16 I 20.0 (5.6), TVimb = 8/6.0 (0.0), AWsan = 5/2.5 (1.0) 
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where it is assumed that <&p M = <3?p r . The expression of the electron neutrino survival probability P, 
that keeps into account the matter effect [17\ 118]. is: 

{U?i for normal mass hierarchy 

(13) 
Ug^Pf + Ue$(l — Pf) for inverted mass hierarchy 

where we have to distinguish the two arrangement of the neutrino mass spectrum compatible with 
present knowledge of neutrino oscillations (see e.g. [H]). We adopt the conventional decomposition of 
the mixing elements in terms of the mixing angles: U e3 = sin #13 and U e \ = cos 9 12 cos #13. We see that 
in the case of normal mass hierarchy, the probability P ~ 0.7 is reliably predicted and rather precisely 
known. Instead, for inverted mass hierarchy, P depends strongly on the unknown mixing angle #13. In 
fact, the so called flip probability Pf (that quantifies the loss of adiabaticity at the 'resonance' related 
to the atmospheric Am 2 ) is: 



P f (E v ,6 13 )=exp 



U 2 3 /20MeV\ 2/3 



3.5 x 10~ 5 V E v J 



(14) 



where the numerical value corresponds to the supernova profile N e ~ 1/r 3 given in [18]. For the 
measured solar oscillation parameters, the Earth matter effect is expected to be pretty small [20], and 
we will neglect it in the rest of the analysis. 

Results for normal hierarchy A fit to the data, made in the same way as for eq. [9] but accounting for 
oscillations with normal mass hierarchy yields: 



M a = 0.5 M , T a = 2.1 ± 0.1 MeV, r a = 0.70tg;^ s 
R c = 13+1 km > T c = 5.1±g:? MeV, t c = 4.4lJ;f s 



(15) 



The correlation coefficients p(x,y) = a(x,y)/(a(x)a(y)), in percent (%), are: 

p(T c , R c ) = -89, p(T c , t c ) = -42, p(T a , r a ) = -16, p(T a ,r c ) = 11, p(r c ,R c ) = 9, 
p(R c ,T a ) = -9, p(r c ,T a )=4:, p(T c ,r a ) = -3, p(T c ,T a ) = 1, p(R c ,T a ) = Q. 

The larger ones correlate the parameters of the cooling phase (a well-known result). The coefficients 
that correlate the phases of cooling and accretion are instead small. 

The x 2 increases by A% 2 = 0.8, namely, it does not change significantly. In other words, we claim 
that in the context of the discussion SN1987A data alone do not provide us with a strong sensitivity 
to oscillations with normal mass hierarchy. This can be understood as follows. The main impact of 
oscillation is on the neutrinos emitted during accretion; in fact, the flux of antineutrinos is multiplied 
by P ~ 0.7. This means that the expected number of events during accretion is the same that we 
found without oscillations, for a value of M a = P x 0.5 ~ 0.35. From table [1] and figure [U we see 
that this case has a perfectly acceptable x 2 5 that concludes the explanation. However, the fraction of 
energy emitted during accretion does not diminish with the oscillations, because even if 30% of the 
antineutrinos become invisible to inverse beta decay reaction, they should be anyway produced. In 
fact, the higher temperature T a implies that f a increases: 

= 8.1 x 10 52 erg . , „ , . 

-- 1.7 x 10 53 erg /"( osc - n - h -) = 32% ( 16 ) 

compare with eq. [TUJ Eq. [10] is the main result of our analysis, and will be discussed in detail later. 
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Figure 2: The x 2 f or the hypothesis of inverted mass hierarchy and for values of #13 allowed by 
present knowledge of oscillations. See text for discussion and warnings. 

Results for inverted hierarchy The case of inverted mass hierarchy can produce P ~ and thus a 
strong suppression of the electron antineutrino flux during accretion. This can be seen clearly from 
figure [21 where (fixing M a = 0.5M & ) we give the value of the x 2 f° r various values of #13 that are 
allowed by what we know at present on oscillations. Thus, the conventional treatment of oscillations in 
eq. [13] implies that the high-luminosity feature visible in the first second of emission disfavors the case 
when neutrinos have an inverted mass hierarchy and values of #13 larger than 0.5° — 1°. It should be 
noted that for the largest value of #13 the fit improves, but the price to pay is a very large amount of 
energy emitted during accretion; e.g., for #13 = 10° we need £ a = 7.7 x 10 53 erg > £ c = 2.1 x 10 53 erg! 
This can be seen as an indication against the case of inverse mass hierarchy and relatively large values 
of #13. This result is surely interesting, but in our opinion it should be taken with caution due to: 

1. Limited astrophysical information. Indeed, in presence of a relevant muon/tau antineutrino emission 
during accretion, the effect of oscillation weakens (at the same time the fraction f a increases, see eq. [2] 
and discussion therein). 

2. Possible modifications of physics of oscillation. Indeed, neutrino contribute to weak ('matter') po- 
tential, especially during accretion, making the problem non-linear |21j . 

Our treatment of oscillation conforms to the conventionally accepted framework and develops in the 
context of eq. [2j This is the best that we can do at present; once the modifications from points 1. and 
2. above will be precisely quantified, it will be a straightforward exercise to repeat these steps to know 
whether these results change significantly. 
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However, it is interesting to try to explore what could happen by relaxing the assumptions we 
followed till now, and this is what we attempt in the last part of this section. For definiteness, let 
us ask what happens if P = k with k = 0.1 or = 0.5. The fit is similar to the one corresponding 
to M a = P x 0.5A/q in tabled] (thus acceptable on the basis of x 2 ) however we would need to emit 
during accretion an energy 1/P times larger than the one we read from table [TJ More in details, we find 
£ a = 3.3 x 10 53 erg for k = 0.1 and £ a = 1.0 x 10 53 erg for k = 0.5. The first case can be questioned on 
theoretical basis, whereas the second, that leads to f a = 36%, seems much more reasonable. This means 
that if non-linear effect are able to produce an effective survival probability P ~ 0.5, inverted neutrino 
oscillations will be not in serious disagreement with SN1987A observations. Comparing with eq. 1121 we 
see that the two case here considered are similar to the case P ~ (i.e., oscillations conforming to the 
conventional expectations) but with a flux of non-electronic neutrinos during accretion = k x $p e . 
However, the estimated values of £ a should be increased by the factor 1 + 2k due to the presence of non- 
electronic species during accretion. We conclude that within the conventional theoretical framework, 
we can have a good fit of SN1987A neutrinos only with important modifications of the hypothesis on 
the flavor composition of the neutrino flux during accretion (eq. [2]) and/or with a strong modification 
of the oscillation probabilities shown in eq. [T3j These effects could save the inverted hierarchy case, 
but would not change the conclusion on f a , that can only increase in comparison with the case when 
oscillations are neglected and presumably also with the case of oscillations with normal mass hierarchy. 

3.3 Expected range of f a and stability of the result 

From here on, we focus on the case of oscillation with normal hierarchy; similar conclusions apply in 
the case of inverted with small #13 < 1° (that becomes indistinguishable from the previous one for very 
small values of #13) and also without oscillations. 

We evaluated the expected range of f a in eq. [16] as follows: 

i) First we found the impact of M a . By considering the % 2 (M a ) allowed by the data (as in fig. [1]) we 
obtain a corresponding range for f a . We find that f a > 18% at 95% C.L. 

ii) Next, we found the effect of the other 5 astrophysical parameters. By propagating the errors, we 
obtain the absolute error 8f a = 10%. This means that f a > 16% at 95% C.L. 

The large errors are due to the fact the the number of events occurring during accretion is only about 100 
We also checked the stability of our results on f a : 

a) By shifting the offset times within their error-bars. 

b) By considering alternative, reasonable parameterizations of the antineutrino flux such as those pro- 
posed and studied by [7]. 

c) By removing some event from the dataset (e.g., when this is done with the event number 10 of 
Kamiokande-II-see [9]-the \ 2 somewhat improves). 

In none of these cases, however, our conclusion on f a does change significantly. The reason is simply 
that this result is not due to an ad hoc statistical analysis, but to two clear experimental facts: 

1) there is a relatively large number of events in the first second, especially in Kamiokande-II dataset; 

2) their energy is relatively low: compare with Sect. 3.2.3 of [9]. 

In conclusion, the data suggest that a considerable fraction energy was emitted during accretion. 
Though it is possible to have a smaller f a for certain values of the astrophysical parameters of the 
neutrino emission that permit a reasonable description of the data, the fact that these parameters are 
not known at present and the occurrence of oscillations suggest that f a was 20% or larger. 

4 This will be certainly much larger when we will observe the neutrinos from a galactic supernova. 
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4 Discussion 



We showed that data from SN1987A imply at 3.2a the existence of an initial phase of intense neutrino 
luminosity, that resembles what is expected during accretion. The agreement of the various astrophys- 
ical parameters with expectations is rather good, but there is a hint that the amount of energy emitted 
during accretion is a few times larger than what is expected in standard calculations. The significance 
of this hint is limited by the small number of detected events, but at the same time, it is a stable 
outcome of the analysis and it is due to clear features of the data (especially, those of Kamiokande-II). 
This result is not contradicted but rather reinforced if three flavor oscillations of neutrinos occurred 
with normal mass hierarchy (see eq. 1161 and discussion therein). Instead, it is not easily reconciled with 
oscillations, if the mass hierarchy is inverted and #13 is larger than about 1°. 

Let us comment these results. The first implication regards the nature of the explosion. We feel 
motivated to propose the speculation that a value f a ~ 20 % (or larger) is a characteristic aspect of 
successful supernova explosions. It will be interesting to see whether such an expectation, driven by 
SN1987A data, will meet the findings of future successful simulations of supernova explosions: see [22] 
for a status report. 

The second implication regards neutrino oscillations. We showed that a large oscillation effect on 
electron antineutrinos (i.e., a small survival probability P in eq. I12j> is disfavored by the data, since it 
would tend to dilute the number of events expected in the first second. This result should be taken with 
caution, since it is possible to argue that the expectations on flavor partition during accretion are not 
completely reliable, and/or that the oscillation could have a non-conventional character [21]. We believe 
that it is urgent to clarify this issue, for such a result could have important implications for several 
future experiments (such as long-baseline experiments, double beta decay search, cosmology j!9j). 

Of course, the most important task is left for the future supernova neutrino experiments, that 
should study precisely the first second of supernova explosion: in particular, the number of the events 
and their energy. It would be interesting to monitor the total flux of neutrinos in the first second, 
e.g. counting neutrons in a heavy water detector as SNO by I'D — > vpn. However, it seems fair to 
conclude that the traditional method of investigation (namely, the observation of supernova electron 
antineutrino by inverse beta decay in scintillators or water Cerenkov detectors, possibly allowing for 
improvements in neutron detection) has still a bright future. 

In conclusion, we presented a state-of-art analysis of the first second of SN1987A neutrino emission. 
We hope that this can be useful to understand better what happened in this crucial moment of neutrino 
astronomy and, possibly, to make further steps toward the theory of supernova explosion. 
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A Peculiarities of the analysis 

As recalled, our analysis is similar to the one of Lamb and Loredo [7], with whom we agree within 
errors when we strictly stick to their procedure. The analysis that we eventually adopted in this paper 
departs from their one for the inclusion of oscillations, and for some technical points that we describe 
here: 
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1) Background: Lamb and Loredo fold the measured background curve with the distribution of the 
energy distribution of the events. However, this has the effect of double counting the detector-dependent 
effects on energy measurement, for which the detected energy is distributed around the true energy 
('smearing'). Thus, we prefer to directly use the background curve, setting Bi = B(Ei); namely, we do 
not perform any folding. This has the consequence that the events of Kamiokande below 7 MeV have 
a higher background rate, and those above 9 MeV a lower background rate, while the other ones stay 
almost unchanged. The changes for Baksan are instead negligible. See [9], Appendix A. 

2) Cross section: When evaluating the signal rate S (eq. [5]) we adopt the inverse beta decay cross 
section calculated in [23]. We use the expression for the cross section da /d cos 9 given in eq. (20) there. 
The energy of antineutrino is given in term of the positron energy E e and the angle 6 between the 
antineutrino and the positron directions: 

Ev = E * + 6 (17) 

1 - (E e -p e cos 0)/m p ' 

where 5 = {jn 2 n — m 2 — m 2 )/(2m p ) 1.294 MeV. We replace cos# in the previous equation with the 
measured values for the first 12 Kamiokande-II and for the 8 1MB events, and set instead cos 6 = for 
the 5 events of Baksan and the last 4 events of Kamiokande-II. 

3) Efficiency: Following [7], we average the signal S as a function of the true value of the energy E e 
over its distribution (assumed to be Gaussian); namely, we keep into account the energy smearing of the 
signal: see eq. [71 But differently from [7] , we include also the detection efficiency as a function of the 
true energy of the event: see eq. El With this procedure, we describe the fact that the expected numbers 
of signal events (i.e., the crucial input of the likelihood) should include all relevant detector dependent 
features (such as loss of events due to light attenuation, fluctuations of the number of photoelectrons, 
detector geometry, etc) including those that lead to an imperfect (% < 1) detection efficiency!! We 
argue in favor of our procedure by considering the following situation. Imagine two detectors where the 
signal interaction rate is equal to the background rate, that differ by the detector efficiency: n = 100% 
in the first one, rj = 10% in the second. Now, suppose that each of them observes an event. According 
to the procedure in [7], the probability that the observed event is due to a signal is 50% in both detectors 
(see table IV and discussion therein), that we find paradoxical. Instead, adopting our procedure the 
probability that the event is due to signal is 50% in the first detector and 9% in the second one, which 
seems to us more reasonable. 

Impact of the modifications: The most important effect is the inclusion of the efficiency of detection, 
followed by the new cross section and finally by our assumption on the background. E.g., adopting the 
procedure of [7J the temperature parameter T c , the radius R c and the time constant r c of an exponential 
cooling model (as in eq. [3] but with M a = 0) are 3.7 MeV, 44 km and 4.4 s. When we include the 
efficiency, they become 4.2 MeV, 30 km and 3.9 s. When we use also the new cross section (eq. (25) 
in [23]) they become 4.6 MeV, 26 km and 3.7 s, and when we include the dependence on the cos# 
(eq. (20) in [24]), these values become 4.5 MeV, 27 km and 3.8 s. Finally, with the new background all 
quantities stay practically unchanged. 



5 We have in mind an 'average efficiency' evaluated by a MC procedure, namely 1) by simulating several events 
with true energy of the positron E e but located in the various positions and emitted in all possible directions, 
and then 2) counting the fraction of times that an event is recorded and finally 3) deducing also the smearing 
on the energy (=average error as a function of E e ). For an even more refined analysis of SN1987A data, one 
should evaluate for any individual event the specific detection efficiency and background rate [7J [9]. In our 
understanding, such a correction on individual basis was performed only to assess the errors on the energies of 
the events, see [4]. 
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